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Abstract 

Volatility dynamics of wavelet - filtered stock price time series is studied. 
Using the universal thresholding method of wavelet filtering and a principle 
of minimal linear autocorrelation of noise component we find that the quan- 
titative characteristics of volatility dynamics of denoised series are noticeably 
different from those of the raw data and the noise. 
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1 Introduction 



Distinguishing valuable information from noise is one of the major problem 
one faces in developing a quantitative picture of price dynamics. Generically 
one divides market participants into those making rational choices based on 
some adopted strategy and those making effectively random decisions (noise 
traders). Understanding the impact of rational strategies requires therefore 
filtering of observed price series with respect to the noise contribution. The 
main difficulty in separating meaningful signal from noise lies in the im- 
possibility of giving a rigorous definition of what noise is. Usual intuition 
visualizes noise as a high frequency uncorrelated dressing on top of a some- 
what slower evolving signal containing useful information. The problem of 
making this picture quantitative is twofold. First, noise can exist not just on 
some single scale, but on multiple scales. In this case one should denoise the 
raw signal at several potentially relevant scales. Second, having noise with 
vanishing linear autocorrelation does not imply that its higher order nonlin- 
ear autocorrelations are zero. A useful instrument for performing multiscale 
denoising is a wavelet - based filtering, see e.g. [HEJE]. As to the existence of 
high order nonlinear correlations in the noise component, this will be a part 
of discussing volatility dynamics in the present paper. Let us also mention 
that recently wavelet - based denoising of financial time series was discussed 
in [HI CD], where a separation in the characteristics of active (turbulent) and 
inactive (laminar) periods in the filtered signal was studied. 

In the present paper we shall concentrate on studying volatity dynamics 
of wavelet-filtered financial series. Volatility (magnitude of price increments) 
is one of the key characteristics of financial time series. In contrast with 
the classical diffusion (Brownian random walk) where volatility is constant, 
financial time series are characterized by an intermittent pattern of high 
and low volatility periods. Such behavior points to a presence of long-range 
memory of the magnitude of price increments (volatility), see e.g. [HE]. In 
particular, long memory property of volatility dynamics shows itself in a very 
slow power-like decay of volatility autocorrelation^. Proper account for the 
impact of varying time horizons on volatility dynamics allows to construct 
a parsimonious model of stock price evolution allowing to describe many 
important properties of observed stock price evolution j6j [8] . 

In the present paper we focus on studying the autocorrelation properties 
of high-frequency volatility. We analyze the data for 5-minute returns of the 
index MICEX10INDEX and five most liquid stocks EESR, RTKM, LKOH, 

3 Detailed studies of volatility correlations have revealed several characteristic time 
scales staying behind apparent long-range memory of volatity dynamics [6l [7] . 
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SBER, SNGS traded at Moscow stock exchange MICEX during the years 
2003-2005. 



2 Volatility dynamics of wavelet - filtered price 
series 

It is well - known that dynamical properties of stock price series are very 
close to those of the random walk. This refers in particular to a nearly absent 
autocorrelation of stock price incremements. A more detailed analysis shows 
however, that zero autocorrelation does not mean that the consequent price 
increments are independent. This dependence shows itself through more 
complex nonlinear correlators, e.g. those of lagged absolute values of price 
increments. 

For all the instruments under consideration we shall study a series of 
normalized 5 - minute logarithmic returns with zero mean and unit standard 
deviation: 



Filtering of the original series {r} corresponds to its separation into fil- 
tered and noise components: 



As mentioned in the Introduction, denosing financial time series is a cru- 
cial step in unraveling the probabilistic dependence patterns hidden in it. 
The main problem is that a rigorous definition of noise based on some fun- 
damental principles of theoretical finance does not exist. Therefore when 
considering a problem of filtering noise from a series an operational definition 
of noise should be formulated. In financial applications the most important 
characteristic of noise is, ideally, an absence of any type of predictability cor- 
responding to complete probabilistic independence of the terms in the series 
under consideration. In practice it is practically impossible to have a denois- 
ing procedure ensuring that all nonlinear autocorrelation functions are zero. 
We shall therefore restrict our consideration to the usual linear autocorrela- 
tion function with lag 1 and require that the lagged autocorrelation of noise 
component of price increments pjv(l) vanishes: 



The property (J3j) ensures that the noise component does not contain linear 
predictability, i.e. that a conditional mean (rjv(i + AT)) rN ( t ) = pn ■ r(t) 





r = 



r F + r N 



(2) 




(3) 
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vanishes for any rjy(^) so that noise does not contain trivial profitable oppor- 
tunities. 

The filtering procedure we employ is based on the universal threshold 
wavelet filtering [TTl [12] . It is realized in three stages: 

• Wavelet transformation of the original time series of price returns r(t) 
with discrete wavelet transformation (DWT): 

r(t) w j:n = 2~ j/2 J dt r(t) ip(2~H - n) (4) 

where tp{t) is a mother wavelet function which we will choose as that 
of the Daubechy 2 wavelet. 

• Separating the set of wavelet coefficients w^ n into the noise w^ n and 
signal w F n contributions through universal thresholding procedure [HI 
[T2] where the signal components satisfy, at each resolution level j, the 
requirement 

| w F n \ > \J2 log (5) 

n is a number of wavelet coefficients at level j and o^ n is their standard 
deviation. 

• Inverse wavelet transformation of the set {w F n } producing the filtered 
time series of returns r F (t) and, by default, a noise component r N (t) = 
r(t) - r F (t). 

As described above, to completely specify the optimal filtering procedure we 
need to give an operational definition of noise. Below we shall define opti- 
mal filtering as that leading to the smallest 1-lag autocorrelation coefficient 
of noise component within the adopted filtering procedure. A natural pa- 
rameter within the universal thresholding scheme one can tune is a number 
of levels involved in thresholding process. The considered time series con- 
sist of 16384 = 2 14 5-minute intervals, so one has 14 levels of resolution. 
A customary choice is applying thresholding to all levels starting from the 
third. We have considered a starting level as a tunable parameter allowing to 
achieve optimal (i.e. leading to smallest one-lag autocorrelation coefficient 
of the noise component) filtering. An analysis showed that for the index 
MICEX10INDEX it was optimeal to start with the 7-th level, whereas for all 
five stocks the optimal starting level is the 3-d. An illustration of the effect 
of filtering is given in Fig. 1 in which the raw cumulative normalized returns 
and their wavelet-filtered counterparts for the index MIXEX10INDEX and 
the stock LKOH are plotted. In Fig. 1 we clearly see a difference that a 
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choice of the starting level makes: filtering of the index series starts from a 
much finer resolution than that of a stock. 

Let us now compare some properties of the raw and filtered time series 
r(t) and r F (t). Two most interesting questions are: 

• Does filtering induce linear autocorrelation in the filtered time series? 

• Is the filtered time series more/less non-gaussian than the original one? 

• What are the correlation patterns characterizing the filtered series as 
compared to those of noise and original noisy series? 

The answer is the first question is given by the lagged autocorrelation func- 
tion p(l). To answer the second one one has to introduce a quantitative 
measure of the distance between the distributions of returns V(r) and V(r F ) 
and their gaussian distributions having the same standard deviation. A con- 
venient quantitative measure of this distance is the (normalized) hypercum- 
ulant [I]: 

For the gaussian distribution k — 1, whereas k < 1 means that the corre- 
sponding distribution is leptokurtic (i.e. having positive anomalous kurtosis 
and heavy tails). The corresponding characteristics for all instruments under 
consideration are shown in Tab lei: 



Table 1 



Instrument 


P(l) 


Mi) 


K 


Kp 


MICEX10IND 


-0.0036 


-0.0512 


0.7872 


0.3178 


EESR 


0.0013 


0.0152 


0.8038 


0.203 


RTKM 


-0.0431 


-0.0648 


0.8533 


0.1773 


LKOH 


-0.0128 


-0.0246 


0.8641 


0.2089 


SNGS 


-0.0602 


-0.2872 


0.8856 


0.19 


SBER 


-0.0755 


-0.2618 


0.835 


0.1747 



We see that the autocorrelations of the original series are somewhat am- 
plified in the filtered ones but remain at the insignificant levels (except for 
SNGS and SBER). A much more pronounced effect is seen in the drastic 
difference between the values of normalized hypercumulant k in hte raw abd 
filtered series. We see that the filtered series are much more non-gaussian 
than the raw ones. 
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To address the third question let us now turn to the analysis of the volatil- 
ity dynamics of raw, wavelet-transformed and noise time series. We will use 
a simple definition of volatility as an absolute value of price returns. A fun- 
damental characteristics of volatility dynamics by is a lagged autocorrelation 
function 

m ((|r B |-(|r|))(|r n+fc |-(|r|))) 

One of the most important properties of financial time series is a slow power- 
like decay of g(k) with k showing that a stochastic process governing temporal 
evolution of volatility is a long-range memory one, see e.g. jl]. 

The volatility autocorrelation functions of the raw, filtered and noise se- 
ries for the index MICEX10INDEX and the corresponding averaged volatility 
autocorrelations of five stocks under consideration are shown in Fig. 2. All 
of them indeed demonstrate a clear powerlike behavor 

g(k) = g k~ a (8) 

The slopes {a} of the powerlike fits to the curves shown in Fig. 2 are given 
in Table 2: 



Table 2 



Instrument 


Raw 


Filtered 


Noise 


MICEX10INDEX 


-0.19 


-0.52 


-0.13 


Stocks 


-0.19 


-0.68 


-0.18 



A striking feature of the plots in Fig. 2 and, correspondingly, the values in 
Table 2 is that the while autocorrelation properties of the raw and noise series 
are practically the same, the filtered series has approximately a two times 
larger amplitude g and a steeper decay with a = 0.67 for stocks and a = 0.52 
for the index compared to the slope of a = 0.19 for both of the original time 
series. This shows that wavelet denoising reveals volatility dynamics very 
different from that seen in the noisy orignal series^. This property can be 
very important in terms of examining real degree of predictability in volatility 
dynamics which is of great interest in many practical applications. This fact 
constitutes a main result of the present paper 

4 Let us mention that the sensitivity of calibrating volatility dynamics to noise effects 
is very large, see e.g. the comment in [S]. 
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3 Conclusions and outlook 



Let us formulate once again the main conclusions of the present study. We 
used a procedure of filtering the original price return series by the universal 
thresholding of the coefficients of its discrete wavelet transform. We found 
that 

• The filtered series is much more non-gaussian as the original one. 

• The filtered series is characterized by a drastically different volatility 
autocorrelation function with larger amplitude and slope than in the 
original series. 

• The volatility autocorrelation function of noise is very close to that of 
the initial series. 

In the present paper we considered only one particular characteristics of 
the wavelet - filtered - its volatility autocorrelation functions. Recent studies 
have revealed a number of striking dependence patterns characterizing high 
frequency dynamics of stock prices [131 El HS1 US] • It will be very interesting 
to see how these dependence patterns change when one considers the denoised 
price series. Work on these issues is currently in progress. 
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Figure 1: Raw versus wavelet-filtered price dynamics for 
MICEX10INDEX; (b) LKOH 
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Figure 2: (a) Volatility atocorrelations for filtered, original and noise se- 
ries for MICEX10INDEX. (b) Average volatility atocorrelations for filtered, 
original and noise series for stocks. Dashed lines: powerlike fits for filtered 
autocorrelation; dotted lines: powerlike fits for raw autocorrelations. 
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